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Abstract 



We present a new realization of asymmetric dark matter in which the dark matter 
and lepton asymmetries are generated simultaneously through two-sector leptogenesis. 
The right-handed neutrinos couple both to the Standard Model and to a hidden sector 
where the dark matter resides. This framework explains the lepton asymmetry, dark 
matter abundance and neutrino masses all at once. In contrast to previous realizations 
of asymmetric dark matter, the model allows for a wide range of dark matter masses, 
from keV to 10 TeV. In particular, very light dark matter can be accommodated with- 
out violating experimental constraints. We discuss several variants of our model that 
highlight interesting phenomenological possibilities. In one, late decays repopulate the 
symmetric dark matter component, providing a new mechanism for generating a large 
annihilation rate at the present epoch and allowing for mixed warm/cold dark matter. 
In a second scenario, dark matter mixes with the active neutrinos, thus presenting a 
distinct method to populate sterile neutrino dark matter through leptogenesis. At late 
times, oscillations and dark matter decays lead to interesting indirect detection signals. 
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1 Introduction 

The neutrino masses, the baryon asymmetry of the universe, and the existence of Dark 
Matter (DM) are the three experimental facts that clearly point to physics beyond the 
Standard Model (SM). Interestingly, it is plausible that all three are related. On the one 
hand, neutrino masses suggest the existence of heavy sterile neutrinos, whose decays in the 
early universe can naturally produce the baryon asymmetry via leptogenesis (for reviews 
with further references see e.g. [21 |3]). On the other hand, the baryon and DM energy 
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densities are of the same order, Q-Du/^h ~ 5, suggesting that they may have a common 
origin. 

One framework that relates the baryon and DM rehc densities is Asymmetric DM (ADM) [U 
|5]. In this framework, the DM particle is distinct from its antiparticle and carries a conserved 
quantum number. An asymmetry in the particle-antiparticle number densities is generated 
in the early universe. Subsequently, the symmetric component is annihilated away by suffi- 
ciently fast CP-conserving interactions, leaving the asymmetric component to dominate the 
relic density. Thus, the relic DM abundance is determined by the asymmetry, rather than 
by the annihilation cross-section, in close analogy to SM baryogenesis, and in stark contrast 
to the thermal DM scenario. 

The ADM scenario has been extensively studied in the literature [HI [TJ El Ej ■ In many 
existing realizations, an asymmetry is first produced in one sector, either in the SM or in the 
DM sector, and is then transferred to the other sector at later times by contact interactions. 
Such a scenario typically predicts similar baryon and DM number densities, ~ nb 

(see however [U [ini [H]), which then allows one to explain the observed Q-Du/^h ratio for 
dark matter mass in the GeV ballpark. However, the lepton and DM asymmetries may be 
produced simultaneously at a very high temperature. This scenario, which we refer to as two- 
sector leptogenesis, may profoundly alter the standard ADM predictions. Several authors 
have previously considered generating the DM asymmetry from leptogenesis [12], however in 
the context of more specific models that predict kdm ~ "^b- In this paper we define a general 
framework for two-sector leptogenesis and demonstrate that it may naturally lead to a large 
hierarchy between ubm and n^. Thus, a wide range of DM masses, from about 1 keV to 10 
TeV, can be obtained within the ADM paradigm. 

The general framework we consider is sketched in Fig. [1} DM resides in a hidden sector 
indirectly connected to the SM via Yukawa interactions with heavy Majorana neutrinos, 
Ni. In this set-up, the SM leptons and the DM particle are charged under an approximate 
lepton number, which is broken by the Majorana masses of N^. The Yukawa couplings can be 
complex, leading to CP violation in the decays of Ni. Throughout this paper we work within 
the framework of thermal leptogenesis. Resonant leptogenesis [13], Dirac leptogenesis [14j , 
or soft leptogenesis [15j could also be considered in this context and would be interesting to 
pursue. 
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Figure 1: A schematic view of our framework: the SM and DM sectors are indirectly connected 
via Yukawa interactions with the same heavy right-handed neutrinos, Ni. The complex couplings, 
Xi and yi, lead to CP violation in Ni decays, and consequently particle- antiparticle asymmetries for 
DM and leptons. 

The generation of the DM abundance adheres to the following steps, 

• A population of (at least) the lightest Majorana neutrino, A^i, is generated in the early 
universe. 

• At temperatures below M^r^ , these neutrinos decay out of equilibrium to both sectors. 
The CP-violating decays lead to a lepton number asymmetry in both the SM and 
hidden sector. 

• As the universe cools well below Mjvj , the washout of lepton asymmetry, and its transfer 
between the 2 sectors, becomes inefficient and the asymmetries are frozen-in. The 
asymptotic asymmetry can, in general, be different in the two sectors due to different 
branching fractions and/or washout effects. 

• As usual, the SM lepton asymmetry is transferred into baryon asymmetry via elec- 
troweak sphalerons. The symmetric baryon component is almost entirely wiped out 
by hadronic annihilations, and only the asymmetric component survives. 

• Similarly, the symmetric component of the DM number density is annihilated away 
in the hidden sector. The relic abundance of DM is set by the remaining asymmetric 
component. DM receives a Dirac mass, m^xX) with another fermion state in the hidden 
sector, X- 
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We present a simple model that realizes the scenario described above. The hidden sector 
contains a chiral fermion x " the DM candidate - and a complex scalar 0. The two are 
coupled to heavy sterile neutrinos via Yukawa interactions, XiNiX(f)- We assume that any 
asymmetry carried by (p is immediately washed out by fast interactions, (p ^ however 
only interactions mediated by the sterile neutrinos can turn x iiito its antiparticle. In this 
way, an asymmetry n^.^ = ~ ^^"^ survive at low temperatures below M^Vi . We also 
assume that additional interactions are present that allow the symmetric component of DM to 
annihilate. For example, there may be a hidden U{l)d gauge symmetry, in which case x ^i-nd 
X can annihilate into hidden photons, X + X — 7d + 7d, as in Ref. [7]. Within this framework 
we study the range of asymmetries, in the SM and DM sectors, that can be generated during 
the 2-sector leptogenesis. In particular, we demonstrate that the asymmetries in the two 
sectors may end up being vastly different, that is to say, ^ 't-al or ^ ^al- In such 
cases, the dark matter mass must be much larger, or much smaller, than a GeV, in order to 
recover I^dm ~ fif,. We argue that, within this framework, DM masses in the keV - 10 TeV 
range are easily obtained without violating any phenomenological constraints. Of particular 
interest is that ADM can accommodate very light dark matter, in the cold as well as in the 
warm regimes. 

The scenario we have outlined may be varied in many ways. In this paper, we discuss two 
simple variations of the hidden sector in more detail. In one realization the cf) asymmetry, 
generated by A^i decays, is not washed out during leptogenesis, but instead survives together 
with the X asymmetry. We will see that the two asymmetries have about the same size, 
nAx — n/\^. Then, at later times, when annihilations are already too slow to significantly 
reduce the DM abundance, decays to x (plus SM states) through an interaction generated 
by integrating out the right-handed neutrinos, and cancels out the DM asymmetry. Thus, 
in this scenario, the DM relic density is set by the primordial asymmetry, but today DM 
is symmetric. This fact has important consequence for phenomenology, notably for indirect 
detection of DM, as DM particles in our Galaxy may annihilate with a large cross-section. 

In another realization, cj) obtains a VEV. Consequently, through Ni interactions, x mixes 
with the left-handed neutrinos of the SM. This scenario is a novel realization of the sterile 
neutrino framework for DM (for a review and references, see [16]), admitting a new and 
simple mechanism to populate its density. The mixing with neutrinos opens up decay modes 
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of DM into SM states. As we show below, for sufficiently small A (0) , x has a long enough 
lifetime to be consistent with current bounds, while the decays may be observable in the near 
future. The VEV of also generates a Majorana mass for x, leading to particle/antiparticle 
oscillations for DM, x ^ X- We show that these oscillations can lead to a large annihilation 
cross-section in our Galaxy. 

The paper is organized as follows. In Sec. |2j we explain the workings of 2-sector lepto- 
genesis. To this end we construct a simple toy model that captures most of the physics and 
discuss washout and transfer effects that influence the flnal asymmetries. Sec. [3] describes the 
complete SM plus hidden sector scenario. In Sec. |4} we discuss the variations of our scenario 
in which, although the DM relic density is set by the x^ X asymmetry, we can repopulate 
symmetric DM at the present epoch. In Sec. |5} we discuss some additional constraints on 
our scenario that arise when the DM particle is lighter than the GeV scale. Our conclud- 
ing discussion appears in Sec. [6j In the Appendix we discuss some technical details of the 
Boltzmann equations of 2-sector leptogenesis. We note that readers who are most interested 
in dark matter model building, and less familiar with the technicalities of leptogenesis, may 



prefer to read sections 2.1 2.3 where we describe the toy model, and the discussion around 



Eqs. 3.2 and 3.3, where we discuss how the generated asymmetries relate to the DM mass. It 
is then possible to skip to sections |4] and [5} where we discuss model building issues. On the 
other hand, readers who are already familiar with leptogenesis may prefer to skip directly 
to section [3| where we discuss the full model. 



2 Toy Model for Two-Sector Thermal Leptogenesis 

In this section we discuss the simultaneous generation of matter-antimatter asymmetries in 
the SM and hidden sector during thermal leptogenesis. In order to highlight the relevant 
physics of 2-sector leptogenesis, we start with a simple toy model. The extension of the 
theory to accommodate the full SM is straightforward, and we highlight the ingredients in 
Sec. El 
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2.1 Toy Model 

Consider two fermion (matter) fields /, x ^i^nd two complex scalars /i, coupled to two Majo- 
rana neutrinos A^j, i = 1,2, 

- £ D ^M,N^ + yMlh + XiNiX<P + h.c. (2.1) 

As the notation suggests, /, /i is a proxy for the SM sector, while Xi 'P represent the DM sector. 
At low energy, the theory admits an approximate global "lepton" symmetry under which / 
and X ci-re charge +1 while Ni are charged —1. The symmetry is exact in the limit Mj — )■ oo. 
On the other hand, we assume here that the quantum numbers carried by h and are rapidly 
washed out by some other interactions. We take both x and to get masses at low energies, 
and DM stability follows from < m^, and an assumed Z2 symmetry under which both 
X and are charged. The DM mass must be Dirac in order to preserve lepton number, so 
we assume that DM gets a mass with another fermion, m^xXi where x has lepton number 
— 1. Furthermore, we assume the presence of additional lepton conserving interactions that 
rapidly thermalize /, h, x, and ultimately annihilate the symmetric component Z + T, x + xQ 
We will discuss these important model building issues later in this paper, but for the moment 
we focus on the mechanism generating asymmetries in the two sectors. The key point of this 
example is to demonstrate that / and x can easily have different asymmetries, allowing for 
a large range of DM masses. 

2.2 Decay Asymmetry 

In order to generate an asjnumetry in leptons and in DM, there must be CP-violation in the 

decays of Ni. While one Yukawa phase in each sector can be removed by field redefinitions, 

the remaining 2 Yukawa phases are physical and lead to CP violation. We use the basis 

where yi and Ai are real and positive while y2 = \y2\G^'^^ and A2 = |A2|e*'^'. We take the 

hierarchal approximation. Mi <^ M2 and assume that we can integrate out A^2 and only 

include A'^i in the Boltzmann equations. We're interested in asymmetries in the decays of 

^The canonical example for such interactions is a U{l)d x Gsm gauge symmetry under which each sector 
is charged separately. Washout of h and (f> asymmetries can be due to Yukawa interactions with other light 
fermions in the theory. 
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r (TVi ^ x0) - r (iVi ^ r (tVi ^ z/^) - r {n, ^ r/^t) 

= r ' ^ r ' ^^-^^ 

I- Ni I- Ni 

where T^^ — {yf + XDMn^/Wtt is its total width. The asymmetries are straightforward to 
compute, 

Ml 1 



^x — 



M2 167r(|/2 + A?) 



(2A2|A2|'sin(2(/)J+yiy2Ai|A2|sin(0i + 0j) , (2.3) 



~ W2 l6Tr{yl + Xl) Wll/2rsin(20O+|/i|/2Ai|A2|sin(0; + 0j) . (2.4) 
We see immediately that the 2 sectors may have different asymmetries, and in particular, 

^ ^ 2r sin(20;) + sin(0; + 0^) ^ ?/i||/2| ^2 5^ 

~ 2r-i sin(20^) + sin(0; + 0^) ' '^~Ai|A2r 

Therefore ej/e^^ ~ r for generic phases. When the couplings of matter fields to both right 
handed neutrinos are similar, yi ~ y2 and Ai ~ A2, the asymmetry for each sector scales as 
the branching ratio of A^i decays into that sector. Of course, when Yukawa couplings within 
one sector are hierarchical, e.g. yi ^ y2 and/or Ai <C A2, the decay asymmetries do not have 
to be correlated with the branching ratios. 

The final asymmetry in each sector is determined not only by the decay asymmetries 
Ex, but also by washout effects and transfer effects that may change the asymmetry in one 
or both of the sectors. These may change the simple dependence on the branching ratios 
quite drastically and may even result with a larger asymmetry in the sector with the smaller 
branching fraction and decay asymmetry. The range of possible asymmetry patterns is 
therefore very rich. 

2.3 Boltzmann Equations 

The cosmological evolution of the sterile neutrinos and the asymmetries are described by 
the Boltzmann Equations (BEs). We introduce the abundance yield — n^/s where 
is the number density of the particle x and s is the entropy density. We are interested in 
the evolution of the asymmetries 1a/,x = ~ ^i,x ^ function of time (or temperature 
T), assuming these asymmetries vanish at early times. To this end we solve the BEs that 
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Figure 2: Feynman diagrams contributing to the 2-to-2 terms in the Boltzmann Eqs. (2.7) and 



(2.8), that transfer the lepton asymmetry between the two sectors. The top row shows diagrams that 



violate lepton number, while the transfer diagrams in the bottom row conserve lepton number. 



include the A^^i decays, inverse decays, and 2-to-2 scattering of matter in both sectors. For 
the initial conditions we assume that the matter in the two sectors is in equilibrium with the 
same temperature, while for A^^i we consider two cases: either equilibrium or zero abundance 
at early times. For the toy model at hand the BEs take the schematic form, 
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(2.6) 
(2.7) 
(2.8) 



Here z = M^i/T, Hi is the Hubble parameter at T = M^^, s is the entropy density, ^iVi,«,x 
are the equilibrium number densities, and Br^ / denote the branching fractions of A^^i into 
the two sectors. Finally, 7_d is the thermally averaged A'^i decay density. 



Id 



m%^Ki{z) 

Tl'^Z 



(2.9) 



Further details and the complete set of equations are given in the Appendix. 
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The first equation describes the evolution of A^i abundance due to decays and inverse 
decays. The strength of these interactions is set by 7^) (or more appropriately T^-^/Hi), 
which thus controls the departure of A^^i from thermal equilibrium. We will always assume 
that Tn-^/Hi is not too small, so that A^^i decays before dominating the energy density of the 
universe. In the other two equations, the terms proportional to ei,^ source the asymmetries 
YaxAi oiice Ni drops out of thermal equilibrium. The terms proportional to Br^^^ describe 
the effect of 2 — )■ 1 inverse decay processes aa ^ Ni {a = l,x or the corresponding scalars) 
which lead to a washout of the asymmetries. The (2^2 washout) stands for AL = 2 
processes aa ^ aa with an off-shell A^^i, while the (2^2 transfer) terms stand for AL = 
or 2 processes, aa <-> bb, bb, that transfer the asymmetries from one sector to the other. The 
Feynman diagrams contributing to these (2-to-2 transfer) terms, are shown in Fig. [2] 



There are two basic regimes of the BEs (2.6 2.8) 



• The narrow- width approximation: Fyvj ^ and T'\j^/ M^^Hi ^ 1. Here, the 
inverse decays are the dominant source of the washout and the 2-to-2 contribution can 
be neglected. Consequently, the last two equations decouple from each other and the 
asymmetries evolve independently in each sector. One can further distinguish the weak 
washout regime, F^v^ ^ M^v^, and the strong washout regime Fat-^ > Mat^. In the latter, 
the asymmetries can be sizably reduced by the inverse decays. 

• The large washout /transfer regime: F^v^ ~ M^-^ or V"]^^/ Mjq^Hi > 1 . Here the 
2-to-2 contributions are important and may change the final asymmetries by many 
orders of magnitude. The last two BEs never decouple, and the 2 asymmetries get 
correlated due to strong washout and transfer effects. 

To describe these washout effects, it is convenient to parametrize the asymptotic asym- 
metries, 

= V.e.Y;,liO) ■ (2.10) 

In the narrow-width approximation, Y"^. is proportional to ex which parametrizes the CP- 
asymmetry of Ni decays, as discussed earlier. In addition, the washout efficiencies, rj^, 
defined through the above parametrization, capture the washout effects occurring during 
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Figure 3: Solutions to the Boltzmann equations for a 2-sector toy model with both sectors in the 
weak washout regime, Tn-^J^ti^^/Hi ^ 1 (and consequently in the narrow-width limit), assuming the 
initial condition IjVi = 1^^^ . In this limit, the washout efficiencies are r]i^y. = 1 and the final lepton 
and DM abundances depend only on ei^^ as in Eq. (2.10). The left plot shows the Ni abundance 



(purple line) as a function of z = M^-^/T, with its equilibrium value, Y^, plotted for reference 
(black dashed). The right plot shows the asymmetry abundances normalized to the asymptotic 
lepton abundance for AL (blue dashed) and Ax with = keV (red dotted) and = 10 TeV (red 
line). 



the thermal evolution. Note however that, while intuitive, these definitions need not imply 



that A^i ever reaches it's equilibrium abundance. In Sec. 2.4 we highlight some effects of the 
2-to-2 scatterings. We comment that in the narrow-width limit, since the BEs for the two 
asymmetries are decoupled, the efficiencies are bounded, rix < 1. However, more generally, 
we will see that 2-to-2 transfer effects can dominate one of the asymmetries, leading to 
rjx > 1. For now, we show in Fig. [3] simple solutions to the BEs for which rji^^ = 1, occurring 
in the narrow- width limit and the weak- washout regime Tn^/Hi <^ 1. The solutions assume 
an initial thermal abundance of A^^i and demonstrates the two viable DM mass limits of keV 
and ~ 10 TeV. 

2.4 Washout Effects 

Washout effects may play an important role in one or both sectors. For a given sector, one can 
distinguish between weak, 'Qi^V jq^/ Hi <^ 1, and strong, Br^Tj^^/Hi ^ 1 washout. We now 
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briefly highlight a few interesting washout effects. To classify the spectrum of possibilities 
it is convenient to divide the discussion into three cases depending on whether each sector 
is in the weak or strong washout regime. 



2.4.1 Weak/ Weak 

In this case Btl^^Tni/Hi ^1. It follows that for M^^ ^ 10^® GeV, the two sectors are in 
the narrow-width regime, Bi l^^Tn-^/Mn-^ <^ 1, and are therefore decoupled. Consequently, 
washout is negligible and the final asymmetry strongly depends on whether A^^i thermalizes 
before decaying. In the case that it does the final asymmetries are set by and one has, 

r/i~l, 77^-1. (2.11) 

This situation is depicted in Fig. [3j The ratio of the asymmetries, 

i?A^iXjr^-, (2.12) 



is then simply Rj\ ~ e^/e^, which can be extracted from Eqs (2.3) and (2.4). In particular, 
for our toy model, Ra — yiy2/^i^2- A hierarchical i?A may appear if the Yukawa couplings 
display a hierarchy between the 2 sectors. 

If Yati (0) = 0, the asymmetry vanishes in the first approximation, as a negative asym- 
metry generated at 2; <^ 1 (when Y^n is still less than ) cancels against the positive 
asymmetry generated at 2; > 1 (when Y^^ > FjVi)- A small asymmetry arises thanks to the 
washout effects being different at small and at large z. One can estimate, 

Vl^^Btl, ^x^§yBr,. (2.13) 

Thus Ra — eiBri/e^Br^, which, again for the toy model, gives yfy2/^i^2- In this case even 
a small hierarchy between the Yukawa couplings of the 2 sectors to heavy neutrinos may be 
amplified into a hierarchical Ra- In Fig. |4]we demonstrate the above scaling with a solution 
in the weak/weak limit for both Yn-^{0) = and ¥^^^{0) = Y^'^(O). 

2.4.2 Strong/Strong 

For large Btl,x^Ni/Hi ^ 1 there is no significant dependence on initial conditions for A^i, 
however there is a qualitative dependence on whether the narrow-width approximation holds. 



11 




Figure 4: Solutions to the Boltzmann equations in the case where both sectors are in the weak 
washout regime Tj\f-^3rL^^/ Hi <^ 1 (which implies the narrow-width approximation), for Br ^ = 10~^ 
and eL,x = 10^^ x Br/,^^. Two solutions are shown, assuming thermal and zero initial conditions for 
Ni. The former implies washout efficiency of order one, r]L,x = 1 as in Fig. On the other hand, 
in the latter case Ni never thermalizes and consequently the efficiencies are smaller, as predicted 



in Eq. (2.13). The left plot shows the ratio of lepton to dark matter abundance as a function of 



z = Mn-^/T while the right plot shows the normalized Ni abundance. 

For Fatj <^ Matj and F^^ ^ Mn^Hi transfer and 2^2 washout effects are suppressed, 
and the asymmetry in each sector evolves independently. One can estimate the washout 
efficiencies as 

Hi I Hi 1 , , 

The asymmetry ratio now scales as Ra — ^l^t^x/^x^^l — ^iy2/yi^2- This is of order one 
unless the couplings to A^^i and N2 display a hierarchy in at least one of the sectors. Thus, 
assuming no hierarchy of Yukawa couplings within each sector in this case, one predicts 
comparable asymmetries in the two sectors. In the ADM scenario that would translate to 
a prediction that the DM mass is comparable to the baryon mass. Thus the strong-strong 
case without Yukawa hierarchies may be linked to the ADM scenarios in the literature that 
predict Yb/YcM ~ 1. In reality, however, the Yukawa couplings in the SM sector typically do 
display a hierarchy to match the observed neutrino masses, which destroys this prediction. 
For F7V1 — M7V1 and/or F^^ ^ M^^Hi washout gets amplified due to the contribution 
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Figure 5: Normalized abundances oflepton and DM asymmetries as a function of z = M^^/T . The 
dashed curves show the expected asymptotic asymmetries for unit washout efficiencies, i]l,x — ^■ 
The left plot shows the solution in the case Tn^/Hi = 1. The DM asymmetry changes sign due 
to the significant washout and transfer effects. The final asymmetry in that sector is greater than 
one. Here Tn^/Mj^i-^ = 0.1 and Br^ = 0.9. On the right plot, the solution is shown with identical 
parameters except for Tj^^/Hi = 50. The corresponding theory is in the strong/strong regime with 
wide Ni width. As can be seen, the large washout and transfer effects reverse the ratio of lepton to 



DM abundance, rendering a larger number density in the dark sector. As discussed in Sec. 2.4 the 
ratio of the asymptotic abundances is independent of ei^^. 



of 2-to-2 processes, which may change the efficiency estimates in Eq. (2.14) by many orders 
of magnitude. See the Appendix for details. The asymmetries in the two sectors can then 



get correlated due to transfer effects and since the transfer terms in Eqs. (2.7), (2.8) are 
proportional to the branching fraction, one expects -Ra oc Brj|,/Bri independently of e^- 
Interestingly, these results imply that far from the narrow-width limit, it is generic to have 
a larger density in the sector with smaller branching fraction. This is demonstrated on the 
right side of Fig |5| -Ra is indeed found to be independent of ex and the density in the DM 
sector dominates. One can also note that the DM density changes sign due to the strong 
transfer effects. 
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2.4.3 Strong/Weak 

An intermediate case occurs when only one sector (say the SM), is in the strong washout 
regime while the other (here the DM sector) is in the weak washout regime. In the narrow- 
width approximation one then finds 

In the first case is suppressed from the naive value cl/^x by a small factor //i/r^ViBri 
(while in the other case the factor is i7^/r|f^BriBr-^ which may or may not be small). 
As in the strong/strong case, the above suppression may allow for a situation where the 
density is larger in the sector with smaller branching fraction. For instance, for eL/e^ — 1, 
Tn^/Hi — 10 and Br^ = 10~^, one finds an order of magnitude larger density in the hidden 
sector. For FjVi <^ MjVi and F^^ ^ M^^Hi the washout of Ya; becomes even larger, further 
strengthening the aforementioned effect, while loosing the ex dependence of the ratio. 



3 Towards a Complete Model 

Let us now briefiy discuss how the above model is modified when we replace one of the 
sectors with the complete SM, 

- £ D + Yi^NiL^H + XiNiX<P + h.c. , (3.1) 

where a — 1 ... 3 counts the SM generation. As in in the toy model, we assume that DM 
pairs up with another fermion, x, to receive a Dirac mass at low energies, m^xX- We define 
the decay asymmetry into the SM as the sum of decay asymmetries into each generation, 
cl — X^Q, eL„. X is now the asymmetric DM candidate. To match observation, the asymptotic 
asymmetries should have the numerical values, 

YrL-eLVLY;,l{0) ~ 2.6x10-1° (3.2) 

yr^-^xVxY^liO) ^ 4xl0-i°(^) (3.3) 

where (0) = 135C(3)/47r^5f* and g^, ~ 100 counts the total number of relativistic degrees 
of freedom at T ~ Mjv^. On the right-hand side, we show the asymmetries that are required 
to reproduce the observed baryon and DM abundances, assuming Ye = 12Yal/37. 
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The value for eiTjL is fixed by experiment (up to a small dependence on (7=,, at T ~ M^vJ. 
On the other hand, a prediction for e^r^^ translates into a prediction for the DM mass required 
to match the observed abundance. For example, if the set-up predicts Y^JY^^ ~ 1, the 
required DM mass is in the GeV ballpark. However, as we discussed in the toy model, such 
a relation between the asymmetries is not a generic prediction of two-sector leptogenesis in 
most of the parameter space: the ratio of the decay asymmetries ^l/^x depends on arbitrary 
Yukawa couplings, and moreover there is a wide spectrum of possible washout efficiencies 
riL,x- All this implies that a large range of dark matter masses are possible. 

If the SM neutrino masses are generated through the see-saw mechanism, a generalized 
Davidson-Ibarra (DI) bound [T7j on M^Vi can be derived. Working in the hierarchical limit, 
Miv^ <^ Mm2^3, "^^^ express eL,x as. 

Mm, Im[(3F*y^ + X*X)M-'YY% . . 

" Stt [2ryt + AA*]n ' ^ ' 

MMjm[{Y*Y^ + X*X)M-^XX*]^, 

" Stt [2yyt + AA*]n ' ^ ' 

where above M = diag(Mjv^, M7V2, M^rg). Concentrating on ei and using the relation 
[yy'"]jj = M]^"^ M]j'^[R iTLi, R^]ij/v'^w with R an arbitrary orthogonal complex matrix, an 
upper bound is found, 

.,<W!e^,„-,CiW. ,3.6) 

Here vew = 174 GeV is the VEV of the SM Higgs and m™**^ is the heaviest neutrino which 
was taken to be 0.05 eV on the RHS. C is a function of M and the Yukawa matrices. It 
is simply expressed in the limit where the A'^i branching fraction into one of the sectors 
dominate^ 

{1 BrL > Br^ 

(3.7) 
iXlMj,JXlMM,y/' BrL<Br^ 

In deriving the small BrL limit above, we assumed is a matrix with order one coefficients, 
to avoid tuning or non-perturbative Yukawa couplings. We see that the standard DI bound 
is recovered when A^^i decays mostly to the SM. A large enough asymmetry in the SM 
requires > 10~^ which then implies Mn, ^ 10^ GeV. In the opposite limit the DI bound is 



^Here Bi'l denotes the total branching ratio of A^i into the sum of the SM flavors. 
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Figure 6: Scatter plots for realistic 2-sector thermal leptogenesis, scanning over models that gener- 



ate the correct SM lepton asymmetry, Eq. (3.2), while at the same time producing the correct spec- 



trum and mixing angles for the active neutrinos. The left plot shows the spread of DM masses as a 
function of the mass of the lightest right-handed neutrino, assuming AIi\f.^/M]\f2 = M^^/Mi^^ = 10. 
The coloring of points indicate the branching fraction of Ni into the hidden sector. 10 keV to 10 
TeV masses are accommodated within the thermal leptogenesis framework. Very light DM is typ- 
ically obtained for smaller M^^ and for hidden sector branching fractions of order 0.1 — 0.9. The 
right plot demonstrates the attainable values for the CP violating parameter, ei, as a function of 
the lightest right-handed neutrino mass. For a given mass, M^i, a maximal value for is clearly 
visible, in accordance with the DI bound. For both plots, we assume that Ni starts with a thermal 
abundance. 



multiplied by a factor depending on the ratio of the Yukawa couplings in the dark sector, and 
on the ratio of the sterile neutrino masses. This typically leads to an additional suppression 
and the bound becomes stronger, unless there is a large hierarchy IA2/A1I ^ M^^/Mn^ 
accompanied by |Ai| > |Yia| > 10~^(M7Vi/10^GeV)^/^ . The modified DI bound is clearly 
visible in the scatter plot in Fig. |6| 

An additional consequence of thermal leptogenesis with see-saw masses is that the SM 
sector typically lies in the strong washout regime, BiiX ]^-^/ Hi ^ 1. Indeed, one finds. 

Hi Svr^t;!^ 0.05 eV ' ^ " ^ 

where ~ 100 and a generic R was assumed (weaker washout may however arise for 
Ri2,Ri3 ^ !)• The hidden sector, on the other hand, may lie either in the weak washout 
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regime for Br^ ^ 1, or in the strong washout regime for Br^^ ~ 1. Note that Br^ ^ 1 
imphes that washout in the hidden sector becomes extremely strong. The muhiphcity of 
available scenarios thus allows for a wide range of the asymmetries that can be generated. In 
particular, it is straightforward to make the dark sector asymmetry subdominant ^ ^zvlj 
corresponding to DM mass larger than GeV. For this the hidden Yukawa couplings have to be 
small, in which case, is suppressed. A theoretical bound on the DM mass, < 10 TeV, 
then arises from the perturbativity bound on the annihilation rate of the symmetric DM 
component]^ 

A lower limit on the DM mass in the 2-sector thermal leptogenesis scenario follows from 
perturbativity, which requires < 10^^. Since Y^^ ~ 4x 10~^, it follows ^ 4x 10"'' and 



Eq. 3.3 yields the lower bound of > keV. Coincidently, the rough astrophysical bound 
on hot DM is also of order keV [19]. In reality, however, keV DM mass is hard to obtain 
in thermal leptogenesis since washout effects typically suppress the initial production of the 
DM asymmetry even if the branching ratio into the hidden sector is large. In principle, such 
washout effects can be suppressed if the two sectors are in the weak-washout regime or if 
the branching ratio into the SM sector is large while at the same time is small. Both 
of these possibilities are harder (but not impossible) to realize in our thermal scenario with 
see-saw neutrino masses, but can easily be found in deformations of this setup. In Fig. |6] 
we show scatter plots that demonstrate the mass reach in the thermal leptogenesis case. In 
these scans we assume hierarchical sterile neutrino masses, M^-^ : M^^ : M^.^ = 1 : 10 : 100, 
which implies < 10~^ and thus a slightly larger lower bound > 10 keV. We see that 
the lower reach of is indeed roughly 10 keV. 

Before closing this subsection, a few remarks are in order. 

• Throughout the paper we have ignored finite-temperature effects which may play a 
significant role in some corners of the parameters space [3j. Nonetheless, we don't 
expect the conclusions to change qualitatively. 

• The conclusions and plots in the above discussion rely strongly on the thermal lep- 
togenesis scenario with neutrino masses arising from the see-saw mechanism. It is 



■^However larger DM masses may be possible when dark matter is composite and there is a hierarchy 
between the confining scale and the DM mass [15] • 
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straightforward to consider other, less hmiting scenarios. For instance the assumed 
hierarchy of the sterile neutrino masses does not need to exist, additional Higgs fields 
may be present, or other leptogenesis scenarios can be the dominant source for the 
asymmetry. In such cases, the DI bound takes a different form and the concluded 
possible DM mass spectrum may be very different. 

• Adding the hidden sector is not exactly the same as adding an additional flavor, as the 
hidden sector comes with a new scalar field. Nevertheless, many effects present in the 
context of 3-fiavor leptogenesis [20] are valid in this case too. 

• In the above we integrated out A^2,3; ignoring the asymmetry produced from their 
decays. As in the SM case, this is not always justified and special care may be needed 
if the Yukawa couplings in one sector display a large hierarchy |21j . 



4 Symmetric Dark Matter from an Asymmetry 

In this section, we consider several simple variations of the framework introduced above. 
In sections |2] and [3} we described the class of models where the decays of a right-handed 
neutrino, A^i, result in a dark matter asymmetry through the operator Nix4>, where is a 
scalar belonging to the dark matter sector, taken to satisfy > m^. We assumed that 
(f) carries no asymmetry due to the presence of fast interactions that convert 0^0. We 
further assumed that (p does not receive a VEV, so that DM is stable. In this section, we 
relax these two assumptions and consider models where: 

1. An asymmetry for is generated. 

2. obtains a VEV. 

We will see below that these simple modifications have important consequences for dark 
matter phenomenology. In models of type (1), where also carries an asymmetry, the decays 
of to X, at low temperature, will reintroduce symmetric dark matter at late times. As we 



discuss in section 4.1 , models of this type predict a large DM annihilation rate in the present 



day. In models of type (2), the VEVs of (0) and the SM higgs, {h), cause DM to inherent a 
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small mixing with neutrinos. In a sense, DM becomes a sterile neutrino with a large Dirac 
mass, my^xX- This has two important phenomenological consequences which we discuss in 



section 4.2 DM can decay into SM fermions giving observable signatures, and x can oscillate 
into X late times, also repopulating symmetric dark matter. 

In both cases discussed here, the symmetric DM component is obtained at late times from 
the dominating asymmetric one. Of course, symmetric DM is also produced directly through 
A^^i decays and may dominate the energy density already at earlier times. The predictions 
of such a scenario are distinct from the ones considered here, and will be presented in future 
work |22]. 



4.1 Restoring Symmetric Dark Matter with Late Decays 

Above, in our example model of thermal leptogenesis, we assumed that (p has interactions 
that are efficient at low energy and set = n^\. This is roughly the situation for the SM 
Higgs, but does not need to hold in the hidden sector. Consider the case where (f) number 
is preserved, so that = — n^t can be nonzero at low temperatures. As we discuss 
below, the asymmetry in results in the restoration of x after x decouples. This leads to 
symmetric DM with some interesting phenomenological possibilities that we discuss in this 
section: a large annihilation rate at the present day and mixed warm/cold dark matter. 

We now outline the cosmology of this scenario. Suppose, as above, that CP violation in 
right-handed neutrino decays produce a x asymmetry, > 0. In the absence of strong 
washout processes, there will also be an asymmetry of 0, n/^^ > 0. In fact, if the BEs for 
the two asymmetries are invariant under the exchange of x and 0, the resulting asymmetries 
are equal: 

riA^ = ■ (4.1) 

One can then check that a sufficient condition for the exchange symmetry to exist, is for 
the low energy interactions to preserve a U{1)l x U{1)^ symmetry under which the fields 
transform as 0(0, 1) and x(l) These symmetries further guarantee that nA^ = 0. 

Incidentally the mass term, m^xx violates the above condition but keeps Eq. (4.1) intact. 



If, however, other interactions exist that involve x, Eq. (4.1) is modified by an order one 



^Here we assume a Maxwell-Boltzmann distribution for the equilibrium values of the two fields. 
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Figure 7: Two-body and three-body decay modes of cf) into DM and SM particles. In the absence of 
washout processes for (j), its asymmetry is equal to that of x- Consequently, the above decays, which 
occur after annihilations have decoupled, repopulate the symmetric component of DM, allowing for 
an observable annihilation signal. 



amount, reflecting the redistribution of tlie x asymmetry into tliat of x- Assuming no such 
interactions, when the temperature reaches the x ci^nd (p masses, the symmetric components 
annihilate and freezeout, leaving asymmetric abundances of x and 0. Then, at temperatures 
below X ciiid decoupling, the scalar can decay to x as shown in Fig. [7} repopulating 



symmetric DM. Eq. 4.1 implies that the final abundances of x and x are equal. 



The Boltzmann equation for Ax now also depends on A0, but given Eq. 4.1, it can be 
expressed in a very similar form to equation Eq. |2.7[ 



—Yax = 71) 



e ( ^Br 



^Ni J ^X 



(2^2 washout + transfer) . (4.2) 



The only difference, from above, is that the 2 — )■ 1 washout term is twice as big as when 
had no asymmetry. 

For simplicity, in the rest of this section, we assume the hierarchical limit Mat^ ^ M^^^.^ 
and only include the lightest right-handed neutrino, A^^i. Furthermore, we suppress lepton 
flavor indices. It is straightforward to extend our discussion to the more general case. In 
our example model, can decay to x through a dimension 5 operator that is generated by 
integrating out the right-handed neutrinos. 

Inserting the Higgs VEV, can two-body decay to x ^i-nd a neutrino. 
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where = ?/^f|^/MAr^ is the see-saw contribution to the neutrino mass matrix, along 
the diagonal in flavor space, from integrating out A^^i. Decays to neutrinos are only weakly 
constrained by Big Bang Nucleosynthesis (BBN), and decays as late as r < 10^ sec are 
allowed if decays entirely through this mode |23J. Decays can easily proceed fast enough, 

_4 /0.1\ V0.05 eV\ / Mjv, \ ^100 GeV 



r,.7xlO-'.ecx(-J ^— -j^^J^— -1. (4.5) 



Equation |4.3| also introduces a three-body decay mode where the Higgs is produced, 
(j) — 7- x^h. This decay is of course suppressed by three-body phase space, but is enhanced, 
relative to the two body decay, when ^ v, 

This three-body decay mode is strongly constrained by BBN because of the hadronic decays 
of the Higgs [211 [25] , implying that (p must decay faster than about a second if the branching 
fraction to three-body decays is appreciable. 

Recall that in the asymmetric DM scenario, x ^"^^ X annihilate with a large enough cross- 
section, aann, such that the symmetric component of DM is subdominant at the present day. 

We must make sure that (p decays late enough such that these x annihilations are decoupled 
(and furthermore, late enough to avoid their recoupling after the decay), so that the x 
abundance is determined by the asymmetry, instead of the annihilation rate. This decoupling 
temperature is defined by the relation, sY^x (o'annv) = H(Tdec), which implies that. 



GeV-^,-/M^V^^), (4. 



- dec 



10 cm^/sec 



100 GeV V i^annV) 

where we've assumed that g^: ~ g^:S at T^ec- Consequently. 



Requiring that (p decays occur after this decoupling temperature leads to a nontrivial con- 
straint on the parameters, 

/ A y / lOO GeV \ / V f lO^ GeV \ ( rn^ 

^* V J\W-''cm^/sJ \ M^, lo.05 eV; ' ^ ' 
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where we've assumed that the three-body decay mode is not dominating, for this estimate. 

If DM is dominantly asymmetric today, then DM annihilations are suppressed at the 
current epoch, making it difficuh to observe indirect signatures of dark matter annihilations. 
But when x is restored by decays of 0, DM annihilation signals are also restored and proceed 
at a boosted rate relative to a thermal WIMP, as in Eq. 4.7[ This provides a novel mechanism 
for producing a large dark matter annihilation rate today. 

In this scenario, such annihilations may account for the leptonic cosmic ray excesses ob- 
served by PAMELA [26], FERMI [27], and HESS [28]. The requirement is simply that DM 
has a TeV scale mass and that its annihilations produce leptons. This provides an attrac- 
tive alternative to the scenario where the DM annihilation rate experiences a Sommerfeld 
enhancement at low velocities p9]. In models with a Sommerfeld enhancement, there can 
be tension between producing the correct relic density and a large enough annihilation rate 
in our Galaxy [30l [31] . But in our framework, this tension is resolved because the DM 
abundance follows from the asymmetry produced by leptogenesis, not from the annihila- 
tion rate at decoupling. We note that both models with Sommerfeld enhancement, and our 
framework, are constrained by the Cosmic Microwave Background (CMB) [32] (and other 
astrophysical constraints, see e.g. [33]) because decays before the time of recombination. 

Finally, we emphasize that decays of provide a mechanism for generating mixed warm / cold 
DM when there is a mass hierarchy: ^ m^. In this regime, x is produced carrying a 
large kinetic energy set by m^. If the decay occurs late enough, such that x is kinetically 
decoupled, then it will not thermalize and there will not be enough time to redshift away its 
kinetic energy. This opens up the possibility for x to constitute cold dark matter while x is 
warm. Parametrically, the velocity of x the present epoch is given by [Mj, 



_5 km f 



s \ 




v^^2^lQ-'—['^\J-^, (4.11) 



where is the lifetime. For example, if = 100 GeV, = 1 MeV, and is one 
second, then x will have a free-streaming velocity of about 1 km/sec, large enough to impact 
the matter power spectrum [31]. These parameters have some tension with the limit of 



Eq. 4.10, but the limit can be satisified at small A and large M^^- This scenario of mixed 
warm/cold dark matter, where half of dark matter is warm and the other half cold, may 
have interesting phenomenological consequences for structure formation, which would be 
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worthwhile to further explore. 



4.2 Asymmetric Sterile Neutrinos from Leptogenesis 

So far, we have assumed that the scalar in the hidden sector, 0, does not obtain a VEV at 
low-energy. This permits DM to be stable when > m^. In this section, we relax this 
assumption by allowing to receive a nonzero VEV. We will see that this simple change 
leads to several new phenomenological possibilities. Because of the nonzero f,^ = (0), DM 
now mixes with the left-handed neutrinos. Therefore, it constitutes a Dirac sterile neutrino 
(we continue to assume that DM has a Dirac mass, which is necessary for its abundance to 
be set by an asymmetry) . This scenario thus provides a novel mechanism to account for the 
correct relic abundance of sterile neutrino DM (for a nice review and references see Ref . [16] ) . 
DM stability is no longer guaranteed, and several decay modes open up due to the mixing 
with neutrinos. For the appropriate DM lifetime, this leads to observable cosmic rays at the 
present epoch. Another consequence of t;^ 7^ is that DM inherits a small Majorana mass, 
fixX^y where fi^ ^ m^. As we discuss below, this leads to oscillations at late times, allowing 



for a large annihilation rate at the present day, as in section 4J^, We consider the above 
effects in detail below. 

Recall that the seesaw Lagrangian (here simplified to the one-flavor case) is given by, 

C D -m^xX + \Mn,N^ + XN^x (0) +yN^L {h) + h.c. , (4.12) 

where have included an explicit Dirac mass for x, and we emphasize that both and the 
SM higgs, h, receive VEVs. After integrating out the heavy right-handed neutrino, A^^i, we 
have the following mass terms, 

CD -m^XX- ~ - f^x'^Xi' + h.c, (4.13) 

where <^ constitutes a small Majorana mass for x, rriu is the usual Majorana mass for 
left-handed the neutrino, and f^^u represents a mass-mixing between x and z/. These masses 
are given by. 
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The X Majorana mass, /x^^, leads to DM particle/antiparticle oscillations x ^ X- Similarly, 
the DM-neutrino mass mixing, /i^j,, leads to DM/neutrino oscillations, x ^- We now 
discuss how these oscillations can modify the cosmological history of this model. 

We begin by considering only the x ^ X oscillations (we will see below that the x ^ 
v oscillations can be neglected for the parameters of interest). It is important that the 
oscillations do not turn on until after DM annihilations decouple, because otherwise the x 
asymmetry, resulting from leptogenesis, is erased. As a consequence, we shall now see that 
the rate of oscillations are slow relative to the Hubble rate at all times. Nonetheless, the 
probability to oscillate becomes sizable at late times, thereby enabling a large annihilation 



rate at the present day, as in Eq. 4.7 This provides an alternative mechanism, compared 



to the late decays discussed in section |4.1[ for generating large cosmic ray fluxes at the 
present epoch [TOt [35]. 

In order to verify that DM does not oscillate too soon, we briefly review the formalism 
for treating particle oscillations in the expanding universe [Ml EZl EH]- Consider a generic 
oscillation of the type a — )■ /3. The BE for production of /3, through oscillations, is given by, 

^ = ^(i^.-/.W)^(l^«-l^/.), (4.15) 

where Pa^p (t) is the probability that a oscillates into /3 after time t, Fq is the total in- 
teraction rate of a, z = m/T is defined in terms of an arbitrary mass scale m, and 
Hi = H{T = m). The oscillation probability, P, is averaged over the interaction time, 
(P) = Va dte~^"^ P. We see that the BE is driven by P x F. Therefore, oscillations 
are in equilibrium whenever PF ^ H and are frozen out whenever PF ^ H . The general 
expressions for P and (P) are: 

Pa^l3 (t) = sm [26^^) sm y-^—tj , {Po^-^P W) = ^ ^^2^ ^ ^2 ' (^'^S) 

where Oap is the mixing angle and Ai^^^ = tosc the energy difference between the states a 
and /3, or equivalently the inverse oscillation time. Note that for simplicity, we neglect the 
effects of finite temperature and density. These corrections are considered, for example, in 
many studies of sterile neutrino production through neutrino/sterile neutrino oscillations |36l 

EZIEH]. 
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We now apply the above formalism to understand x X oscillations. The mixing angle 
and the energy splitting in this case are given by 

0,,^-, Ai?,,^<^ ^ (4-17) 

In order to insure that these oscillations are slow before DM annihilations decouple, it is 
sufficient to require that TP < H aX the temperature of decoupling, 

(P^_,^(t))r^(Tdec) < i/(rdee). (4.18) 

Since the asymmetric component of x cannot be neglected, T^^c must be calculated by com- 
paring the total number density times the annihilation rate to Hubble, Fann = [^asym(^dec) + 
'T-sym(Tdec)] (fann^^) = H (Tdcc) ■ For our casc, whcrc the asymmetric component dominates at 
low temperatures one finds, m^/Tdec ^ 20 implying that even at temperatures below that 
at which the symmetric component decouples, rn^/T^l^ ~ 20, oscillations may recouple the 



annihilation process, thereby significantly altering the DM abundance. Using Eqs. (4.16) 



and (4.17), Eq. (4.18) reads at T = Tdec 



fix<\/^X^-nn{T = T^ec). (4.19) 

This condition can be converted into a limit on the size of Xv^, which enters the oscillation 



probability through /i^. Using Eqs. 4.14 one finds the limit 



^ < 3 X 10- (^X" f 1~Z5) (p^Y' hll I . (4.20) 



^ \10^^ GeV J \ {aannV) J \{(^annV) J ^9*3 

where (Tann is the x + X annihilation cross-section, atot is the overall x interaction cross- 
section (which determines the total DM interaction rate F^), and and g^:S are evaluated 
at the temperature that DM annihilations decouple. We see that the requirement that DM 
particle/antiparticle oscillations are slow, before DM decouples, presents a rather stringent 
limit on the quantity Xv^ relative to the DM mass, m^. 



We note that if Eqs. 4.18 and 4.20, are satisfied, then PT <^ H, also at temperatures 



below DM decoupling, due to the drop in F relative to H as the Universe cools. This means 



that, given the condition Eq. 4.18 is satisfied, oscillations always remain slow in the expanding 
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Figure 8: Dark matter decay modes in the case where {(p) ^ 0. These modes occur due to mixing 
of DM with SM neutrinos and place stringent constraints on the mixing angle for ^0.1 GeV. 
For lighter DM, oscillations into x <ire more constraining. 



universe, and never produce a significant yield of Y^,. Still, there is a large annihilation rate 
today, as long as tosc = f^^^ is much shorter than the age of the universe. In this regime, 
whenever two DM particles collide at the present day, there is an 0{1) probability that one 
has oscillated. We also comment that if fi^ is chosen so that tosc is longer than the timescale 
for recombination, then the constraints from the CMB [32] are alleviated. 

Now we consider oscillations between DM and neutrinos, x ^ Oscillations of this type 
are highly suppressed by the small mixing angle between DM and neutrinos, 

V-*^. AE,^^l'<'^ ^--'"^ (4.21) 



In particular, combining the constraint of Eq. |4.20| with Eqs. |4.14[ |4.16[ and |4.21] we learn 
that, 

This oscillation probability is small enough that PT ^ H for all temperatures. This shows, a 
posteriori., that we were justified to neglect DM/neutrino oscillations in the above discussion. 

There is, however, one very important effect of DM/neutrino mixing, 9^^^. This mixing 
opens up several decay modes for DM. For example, suppose that x is lighter than the 
electroweak gauge bosons. Then as shown in Fig. |8j DM can three-body decay through 
off-shell Z*/W\ X ^{Z* ff) and x ^ e-{W+* -> //'), where ///' are SM fermions. 
The rate for these decays is given, parametrically, by the expression, 

(4.23) 
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Of course, these decays are constrained to not proceed too quickly. In order to be consistent 
with cosmic ray and diffuse gamma observations, > 10^^ s [Ml EH]- This hmit can also be 
expressed as a constraint on the size of relative to m^, 

hi* < 10- (^) ( ^'-^ ] r^^ii^l fl^] (4 24) 

We see that the DM lifetime is a stronger constraint than DM particle/antiparticle oscil- 
lations, when ^0.1 GeV, while oscillations present the dominant constraint on light 
DM. For DM above the GeV scale, the cosmic rays are observable, by ongoing and future 
observations, when Xv^f, saturates the above limit. We note that sterile neutrino DM of this 
type can also decay, at one loop, to a photon line, x '^1 [IB]- However, this decay is 
suppressed, relative to the above three-body decays, by a factor of ~ a /An. We also note 
that in the regime where > m^, DM will dominantly decay two-body to longitudinal 
electroweak gauge bosons and the Higgs, x ^ Zu, x — ^ W^e~ , and x ~^ hu. In this regime, 
the decays proceed faster than above, and the constraint on A is significantly stronger. 



5 Cosmology and Light Dark Matter 

As we have seen, asymmetric DM populated by leptogenesis can accommodate a wide range 
of DM masses, between about keV and 10 TeV. Within this broad framework, cosmology, 
astrophysics, and colliders can all constrain the properties of DM, and the hidden sector 
within which it resides. In particular, the hidden sector must contain additional degrees of 
freedom, lighter than DM, which the symmetric DM component can annihilate into. While 
the constraints on weak-scale DM are well-known (for a review see [30]), light (keV to 10 
GeV) DM raises several interesting issues, which are the subject of this section. 

For concreteness, we specialize to models with a hidden U{l)d gauge symmetry [7j, where 
the symmetric DM component annihilates predominantly into a pair of hidden sector pho- 
tons, 7rf. The hidden sector couples to the visible sector through the vector portal, 

(5.1) 

where F' {F'^^) is the 7^ (electromagnetic) field strength. The above coupling allows 7^ to 
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decay into kinematically available SM fields with non-zero electric charge. The width and 
lifetime of 7^ is determined by its mass m^^ and by the mixing parameter e. 

We will find it helpful to consider two cases separately, depending on whether is 
heavier than, or lighter than, the MeV scale. First, we consider > MeV. In this regime, 
we are also free to assume that m^^ > MeV, and in particular that no sub-MeV hidden states 
exist. This means that there are no new relativistic DOF present during BBN, and the only 
constraint from BBN is that the the lightest state in the hidden sector decays to the SM fast 
enough to avoid late dissociation processes [25] . This requirement is not hard to fulfill. 
For instance, if the hidden photon is the lightest hidden state, it decays to electron pairs 
before T ~ MeV for jH] 

.>3xlO-f^V'\ (5.2) 

Further constraints on '-f^ are summarized in [12]. For a hidden photon heavier than ~ GeV, 
the constraints come from B-factories and muon anomalous magnetic moment and imply 
e < 10^^. For m^^ between MeV and GeV, beam dump experiments (typically assuming a 
2-body decay of the lightest state) require e < 10~^. Finally, for m^^ between ~ 1 — 50 MeV, 
constraints from supernovae cooling imply e < 4 x 10~^. 

Now we consider very light DM, G [keV , MeV]. In this regime, there are two impor- 
tant changes to the above discussion: (1) there are additional light DOF present during BBN, 
and (2) for m^^ < 2me, the hidden photon becomes cosmologically long-lived because there 
remain no lighter states with electric charge. Both of these facts are potentially hazardous, 
and we now discuss how to evade the danger. 

The model can be ruled out by BBN, because DM plus the hidden photon exceed the 
number of relativistic DOF allowed during BBN at the SM temperature. This constraint 
can be avoided if the hidden sector is cooler than the SM when T ~ MeV. This is only 
possible if kinetic equilibrium is not maintained between the two sectors, which requires 
that the following type of reaction be inefficient: 7^6^ ^-t- 76^. This reaction is decoupled 
when T > (e^a;|^^/7r^5fy^)Mpi. Therefore, the two sectors are not in kinetic equilibrium at 
Tbbn ~ 1 MeV if 

e < CBBN ^ 7 X 10-^ (5.3) 
Suppose that this condition is satisfied, and that the two sectors begin at the same tem- 
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Figure 9: The BBN limit on the effective number of degrees of freedom in the hidden sector, g^, 
as a function of the number of degrees of freedom in the visible sector, g^^^ , at the temperature the 
two sectors decouple, T^ec- For reference, the vertical dashed lines indicate the sizes of the full SM 
and MSSM. Here we have assumed that the entire hidden sector is lighter than the BBN scale ~ 
MeV. 

perature (for example at the leptogenesis scale) and decouple at a lower temperature, 
Mat^ > Tdec > Tbbn- If more DOF freezeout in the SM than the hidden sector, before 
BBN, then by the separate conservation of entropy in the two sectors, < T. The limit 
from BBN becomes, at 95% CL 03], 

9'in^^)(P^) <2.52. (5.4) 
V-'bbn/ 

Here T^bi<s is the temperature in the hidden sector at BBN, while g^{T) is the effective 
number of relativistic degrees of freedom in the hidden sector, as defined in |43j. In Fig. [9] we 
use this relation to show the constraint on the size of the hidden sector, g^, as a. function of 
the number of DOF in the SM when the two sectors decouple, (yf^** (T^ec), assuming that all 
hidden particles are lighter than the temperature of BBN. We see that rather large hidden 
sectors can be accommodated, (7^ ~ 50 if the two sectors decouple above the electroweak 
scale. 

The bottom line of the above discussion is that the two sectors cannot be in thermal 
equilibrium at the time of BBN, which enforces the constraint e < esBN- We note that even 
stronger constraints arise from bounds on the lifetime of the sun and horizontal branch stars 
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nH'{T^) = Ch'U'^] n,{T), (5.5) 



that require e < 10~^^ 

Finally, we discuss how to avoid overdosing the Universe, since the hidden photon is 
cosmological long-lived below the MeV scale. One possibility is to drop the assumption that 
the hidden sector has only one mass scale. Then we can make the lightest state in the hidden 
sector sufficiently light, such that it does not overdose the universe even if it is cosmo logically 
stabl^ An estimate for the upper bound on the mass of the lightest state, h' , can found by 
assuming it follows a thermal distribution. At temperatures above its mass, > rrih', the 
number density is related to the photon number density through, 

2 \f 

with Ch' = 3/4 {Ch' = 1) for a fermion (boson) h'. Below m/j/, the energy density scales as 
and is therefore proportional to n^(To). Comparing to the measured DM energy density, 
one finds, 

m,, <3evQ^C,-V, (5-6) 

where ^h' = T^/T at the time when h' becomes non-relativistic. This result is in agreement 
with the more precise calculation of [H]. 

Another intriguing possibility exists. Above we assumed that the hidden sector was in 
thermal equilibrium with the SM at the time of leptogenesis. However if the hidden sector 
couples only weakly to the right-handed neutrino, the two sectors may have never been in 
equilibrium with each other. In that case the hidden sector can easily be much colder than 



the SM, so as to allow for a heavier nih', cf. Eq. (5.6). Furthermore, it is interesting to note 
that if the DM annihilation rate is slow (or even vanishing), the symmetric component can 
dominate over the asymmetric one. Then the symmetric component would be responsible 
for the observed relic abundance, while at the same time it would be related to generation 
of the baryon asymmetry. This possibility is however outside the ADM paradigm, therefore 
we postpone its study to future work [22] . 

As a final remark, we note that the above discussion assumed a hidden gauge group 
coupled to the visible sector through gauge kinetic mixing. Other portals, such as the Higgs 
portal may be considered. The Higgs portal can have very different (and potentially weaker) 



^The lightest state may reside in a distinct sector thus explaining the hierarchy of scales. 
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astrophysical constraints than the vector portal. This is because 7^ couples to electric charge 
and therefore equally to electrons and protons, whereas a hidden sector scalar, that couples 
through the Higgs portal, will couple more strongly to protons than electrons. 

6 Outlook 

In this paper we discussed the ADM scenario in the context of 2-sector leptogenesis. The 
asymmetries in the SM and DM sectors depend on several factors, such as the branching 
fractions and decay asymmetries of the sterile neutrinos, and the strength of washout effects 
in each sector. Consequently, the ratio of dark matter number density to the baryon number 
density is very sensitive to the model parameters, thereby accommodating a wide range 
of dark matter masses from keV to 10 TeV. Our findings suggest that the spectrum of 
predictions of ADM, especially of those concerning the dark matter mass, is much wider 
than previously thought. This is very important for planning direct and indirect experimental 
searches that target the ADM scenario. 

Here, we have chosen to focus on the concrete scenario of thermal leptogenesis with 
hierarchical neutrino masses. There are several variations and open questions that remain 
unexplored: 

• Other leptogenesis scenarios, such as soft leptogenesis, resonant leptogenesis or Dirac 
leptogenesis, may be accommodated within the 2-sector leptogenesis framework. It 
would be interesting to study the above variations as they are expected to admit 
different dynamics and in some cases produce very different phenomenology. 

• It would be interesting to extend our treatment of the BEs to include finite temperature 
corrections and the full flavor structure of the theory. These corrections are well- 
studied for traditional leptogenesis and may have interesting consequences for 2-sector 
leptogenesis as well. 

• In one of the variations studied here, the DM mixes with the active neutrinos, providing 
a novel realization of sterile neutrino DM. It is worthwhile to investigate this scenario 
further. In particular, it would be interesting to understand how the constraints derived 
from late time oscillations or decays can be ameliorated in other leptogenesis scenarios. 
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• It is possible to imagine a similar mechanism for populating ADM in a single-sector 
leptogenesis model [iS]. Clearly the phenomenology of such a scenario would be dis- 
tinct. 

• Another deformation studied here predicts dark matter to be a mixture of cold and 
warm components with roughly equal numbers. This may have interesting phenomeno- 
logical consequences for structure formation, which would be worthwhile to explore. 

• It may also be interesting to consider variants of this scenario where the symmetric 
DM component dominates the relic density. In particular, the symmetric component 
originally produced by the decays of the heavy sterile neutrinos may not annihilate 
away. In such a case, the asymmetric lepton density is suppressed (due to the usual 
bounds on the decay asymmetry into the SM sector) with respect to the DM number 
density, thereby predicting light DM. We postpone a study of the details of such 
scenarios to future work [22] . 
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A Appendix: Boltzmann Equations for the 2-Component 
Toy Model 

For completeness, we discuss the BEs in the 2-component toy model, taking into account 
the 2-to-2 washout and transfer terms. The equations for the asymmetries Iao = Ya — Yd of 
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the components / and x read 



dz 



zKi{z) 
^ Koiz) 



{Y'^\-YN,) + 2Biilw{z)Y^a 



(A.l) 



where a = 6 = x,/, and Hi = H{T = M^rJ. The source term proportional to the 
decay asymmetry of A^^i is the same as in the narrow-width approximation we discussed 
earher. The second term is due to washout processes aa ^ dd within one component, and 
the terms in the second hne are due to transfer processes aa ^ bb and aa ^ bb between the 
2 components. Exphcitly, the washout and transfer terms take the form 



liiz) 



TT 



dteKi{t)f,{e/z^ 



(A.2) 



where f = V]^^/Mjq^ and 

s/2 



fw{s) 



1)2 + r2 



+ 



s - log(s + 1^ 



2(s 



[s + l)log(s + l)-s 



1)2 + P 



, s/2 ^ \og{s + 1) 



s + 1 

s/2 
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Y2 



s + 2 



+ 



s s — log(s + 1) 

I 



s + 1 



3^2 
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Y2 



(g + l)log(s + l)-s (s + 2)log(s + l) -2g 
s + 1 s 



(A.3) 



In the hmit F — )■ the pole terms in fi dominates the integral at small and moderate z, and 
we find Jj ~ z^Ki{z)/4. The BEs then reduce to 

dY^a 



dz 



Hi 



/ N V^nI - ^AfJ + Bra 7. ^Aa 



(A.4) 



K2{z) ^'"1 ' 2 

where the second term describes washout out due to 2 — )■ 1 inverse decays. We see that in the 
limit f — 7- the two components decouple and evolve independently. Their final asymmetries 
are determined by the decay asymmetries and the strength of washout effects in each 
sector set by BiaJ^ n^/ Hi. In addition, the asymptotic value for the asymmetry depends 
on the initial conditions for the sterile neutrino yield Yatj, when washout is weak. See the 
discussion in Sec. 12.41 
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As the pole contributions to the integrals decay exponentially with the asymptotic 
behavior of the BEs is determined by 0{r) off-pole contributions who are only power sup- 
pressed at large z. In other words, at late times 2 — )■ 1 inverse decays become subdominant 
with respect to 2-to-2 off-pole processes. The crossover typically happens for moderate z, 
around Zc ~ 25, where Zc only logarithmically depends on F. For z > ZcWe can approximate 



'^2 r 



W Hi 



ItAz) 
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W Hi 



Ni 



z^ F 



(A.5) 



where W = (32f /7r)(F7Vi/-^i) is the order parameter that sets the strength of the late time 
washout and transfer effects. For z > Zn the BEs reduce to 



dz 



W 



M ■ 



Yai 


M =^ 


Yax _ 





6Br^ + BrjBr^ 



Br^Br. 



Br/Br^ 
mil + Br/Br^ 



(A.6) 



Thus, the asymptotic asymmetry can we written as 

W 



exp 



-M 



Yai{zc) 



Y, 



Ax I 



(A.7) 



in terms of the "boundary conditions" at the crossover scale Zc- It is illuminating to study 
the limit where the branching fraction into one sector, say BR;, is very small, in which case 
the exponent simplifies to 







exp i-—M\ 
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M I ^-6W/zc 
36 



. (A.8) 



For W/zc ^ 1 the diagonal elements are approximately 1, and then late washout does not 
affect the asymmetries. In that case the washout effects are very well described by the 2 — )■ 1 



processes encoded in Eq. |A.4[ The magnitude of the off-diagonal terms describing transfer is 
ByiW/zc ^ 1; nevertheless, transfer could still be relevant in this limit if Yai{zc) ^ ^AxI^c), 
or the other way around. Late washout and transfer become important when W/Zc ^ 0.1, 
that is when TTn^/Hi > 0.1. For Br; ^ 1 we typically expect Yai{zc) ^ Yax{zc) (unless 
ei <^ e-^) because of the 2—7-1 washout at moderate z that is larger in the x sector. If this 
is the case, one effect is the suppression of the asymmetries by the factor e~^'''^/^=, which 
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may or may not be small depending on Br/fTjVi/i?i. Another effect is that the asymmetries 
become correlated, 

ys/yr^ ~ -VBr,, ^ ^ o.i . (a.9) 

Note that the asymmetries end up with opposite signs: particle domination in one sector 
implies anti-particle domination in the other. Amusingly, the larger magnitude of the asym- 
metry survives in the sector with a smaller branching fraction for Ni decays. Contrary 
to naive expectations, transfers do not lead to Y^JY^^ ~ 1, unless Br/ ~ Br^. Thus, 

y^/y^x ~ ^' often advertised as the prediction of ADM, is not a generic prediction of our 
2-sector leptogenesis scenario. 
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